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ABSTRACT 

Photoionization by the high-redshift ultraviolet radiation background heats low 
density gas before it falls into dark matter potential wells, and it eliminates the 
neutral hydrogen and singly ionized helium that dominate cooling of primordial gas 
at temperatures of 10'* — lO^K. We investigate the influence of photoionization on 
galaxy formation using high-resolution simulations with a 1-dimensional, spherically 
symmetric, Lagrangian hydrodynamics/gravity code. We find that the presence of a 
photoionizing background suppresses the formation of galaxies with circular velocities 
Voire ^ 30kms~* and substantially reduces the mass of cooled baryons in systems with 
circular velocities up to Vdrc ~ 50kms~*. Above Vdrc ~ 75kms~*, photoionization 
has no significant effect. Photoionization exerts its infiuence primarily by heating gas 
before collapse; the elimination of line cooling processes is less important. We discuss 
the implications of these results for hierarchical theories of galaxy formation. 

Subject headings: Galaxies: formation. Hydrodynamics, Methods: numerical 

1. Introduction 

Contemporary cosmological theories imply that the galaxies and the large-scale structure 
observed today have evolved by gravitational instability from small-amplitude primordial density 
fiuctuations in the early universe. In theories where the universe is dominated by cold, coUisionless 
dark matter, there is no natural cutoff in the primordial power spectrum at small scales. The 
clustering of dark matter in such theories is hierarchical: low mass objects collapse first, then 
merge into progressively larger systems. However, the material that we observe directly in galaxies 
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is baryonic, so its evolution is influenced by hydrodynamic and thermodynamic processes such as 
shocks, radiative cooling, and photoionization, and these processes can introduce natural scales 
into the baryon distribution. In particular, pressure gradients can suppress gravitational collapse 
on scales below some characteristic mass, usually referred to as the Jeans mass. 

In this paper we study galaxy formation in the presence of photoionization by an ultraviolet 
(UV) radiation background, using simulations of spherically symmetric collapses of mixed 
baryon/dark matter perturbations. The potential importance of photoionization was emphasized 
by Efstathiou (1992; see also Shapiro, Giroux & Babul 1994), who suggested it as a mechanism 
to suppress the formation of dwarf galaxies in hierarchical models. According to semi-analytic 
calculations, hierarchical models naturally predict a faint-end slope of the galaxy luminosity 
function that is much steeper than most observational estimates (e.g. White & Frenk 1991; 
Kauffman, White & Guideroni 1993; Cole et al. 1994). Various solutions to this problem have 
been proposed, including feedback from star formation (Dekel and Silk 1986), reheating of the 
intergalactic medium by supernova winds (Tegmark, Silk, & Evrard 1993), Compton heating from 
energetic objects at very high redshifts (Blanchard, Valls-Gabaud, & Mamon 1992), and rapid 
merging of dwarf galaxies between their formation and z = Q. However, photoionization would be 
an attractively conservative answer to the faint galaxy conundrum because the spectra of quasars 
provide direct evidence that the intergalactic medium was highly photoionized out to redshifts 
z r^b (e.g., Schneider, Schmidt, k Gunn 1989; Webb et al. 1992; Meiksin k Madau 1993). Even 
though some gas must have collapsed and cooled early to produce stars and/or AGNs, the sources 
of the background radiation, it is plausible that most of the gas was still relatively homogeneous 
at the time that it was reionized. Photoionization can have two distinct effects on pre-galactic or 
proto-galactic gas: it heats gas which would otherwise be at very low temperature, and it reduces 
cooling rates by lowering the fraction of neutral atoms. As Efstathiou (1992) pointed out, both 
effects can influence the ability of gas to cool and condense in low mass systems. However, because 
of the simplifled nature of Efstathiou's analytic model, his estimate of the quantitative impact of 
photoionization was very uncertain. 

Thoul and Weinberg (1995; hereafter Paper I) developed a 1-dimensional, spherically 
symmetric, Lagrangian hydrodynamics/gravity code to study the effects of radiative cooling on 
the high-mass end of the luminosity function. In this paper, we use this code to study the effects of 
radiative cooling and photoionization on the formation of low mass systems. While 1-d simulations 
require a geometrical idealization, they complement more computationally intensive calculations 
with 3-d N-body/hydro codes (e.g. Katz, Hernquist k Weinberg 1992; Cen k Ostriker 1993; 
Evrard, Summers k Davis 1994; Navarro k White 1994; Steinmetz k MuUer 1994) in several 
important ways. First, 1-d simulations can achieve much higher spatial resolution and mass 
resolution, so they compute hydrodynamic and thermodynamic effects more accurately. Second, 
with faster simulations one can investigate a broad parameter space, and one can gain physical 
insight by examining how variations of inputs change the flnal results. Finally, 1-d simulations 
are easier to visualize and interpret than their 3-d cousins, again making it easier to gain physical 
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understanding of their results. 

In §2 below, we describe our treatment of photoionization and radiative cooling and discuss 
the initial conditions used for the simulations in this paper. Section 3 presents the results of the 
simulations. We discuss their implications in §4. 



Our simulations evolve a spherically symmetric density enhancement in an Einstein-de Sitter 
(J7 = 1) universe. We assume that coUisionless dark matter and baryons are initially well mixed, 
with a baryon fraction ^If, = 0.1. The coUisionless component is described by the equation of 
motion dv^/dt = —GM(rd)/r'^, where and are the radii and velocities of the coUisionless 
mass shells and M(rd) is the total (baryonic and dark matter) mass inside radius r^. The baryon 
component is described by the fluid equations for a perfect gas moving in the gravitational field 
of both components. We use the standard, second-order accurate, Lagrangian finite-difference 
scheme (Bowers & Wilson 1991). Paper I describes the code in detail and presents many tests on 
problems with analytically known solutions. 



The radiative cooling is calculated for a gas of primordial composition, 76% hydrogen and 
24% helium by mass. We compute the abundances of ionic species as a function of density 
and temperature by assuming that the gas is in ionization equilibrium with a spatially uniform 
background of UV radiation. In other words, we choose the abundances so that the rate at which 
each species is depopulated by photoionization, coUisional ionization, or recombination to a less 
ionized state is equal to the rate at which it is populated by recombination from a more ionized 
state or by photoionization or coUisional ionization of a less ionized state. Ionization equilibrium 
should be an excellent approximation under the physical conditions that are relevant for these 
simulations (see, for example, the discussion of Vedel et al. 1994). The full set of equations that 
we use to obtain abundances and cooling rates is listed in Katz, Weinberg & Hernquist (1995). 
The original source for most of these formulae is Black (1981), and we adopt the high-temperature 
corrections of Cen (1992). 

We model the UV background spectrum as a power-law with a redshift-dependent 
normalization. 



2. 



Model 



2.1. Cooling and heating 




(1) 



where z//, is the Lyman limit frequency (hi'L = 13.6 eV), a is the spectral slope, and 



Jo{z) = J-2i{z) X 10"^^ erg s"^ cm'^ sr"^ Hz 
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Fig. 1. — Cooling and lieating rates as a function of temperature for a gas of primordial composition 
(ffj = nnmp/pg = 0.76, fue = 0.24) and a UV background intensity J_2i = 1 (thick lines) and J_2i = 
(thin lines). The thick dashed lines are the cooling rates, the thick dotted lines are the heating rates, and the 
thick solid lines are the net (|r — A|) cooling/heating rates, all scaled by n'^'^ . Abundances of ionic species 
are computed assuming ionization equilibrium. The units of A and F are ergs s~^ cm~^. 
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is the amplitude. Predictions of the UV background from QSOs and observations of the proximity 
effect in the Lya forest suggest that J-21 ~ 0.3 — f at 2; ~ 2 — 4. In Figure I we show the cooling 
and heating rates F and A and the net cooling/heating rate |F — A| as a function of temperature. 
In the presence of an ionizing background, the cooling curves are also function of density and of 
the intensity and slope of the UV spectrum. In panels (a) and (b) we show the cooling curves 
for uh = 10"'* cm~"^ and nn = 10~^cm~"^, a = 1, and J-21 = 1- In panels (c) and (d), we show 
the cooling curves for the same values of nn and J-21 but with a = 5. Here nn is the hydrogen 
number density, which can be conveniently expressed in the form 



where z is the redshift, pj pu is the ratio of the local density to the critical density pn = (GTrGi^)"*, 
flf, is the baryon density parameter, and h = i7o/(100 km s~* Mpc~*). For comparison. Figure 1 
also shows the cooling curve obtained in the absence of a UV background. 

As Figure 1 illustrates, the UV radiation affects the thermal state of protogalactic gas in 
two ways. The residual energy of photoelectrons heats low-temperature gas, and photoionization 
eliminates the neutral hydrogen and singly ionized helium that dominate cooling at temperatures 
T ~ 10"* — lO^K. The zero-crossings in Figure 1 indicate the equilibrium temperature Teg, where 
cooling and heating balance. Values of Teg lie in the range lO"* — lO^K, depending on the gas density 
and on the UV background spectrum. A harder spectrum yields more energetic photoelectrons 
and a higher Teg. Conversely, raising the gas density with a fixed UV background increases cooling 
and lowers Teg. In the low density limit (tih ^ 10~^ cm~"^ for J-21 ~ 1), Teg becomes independent 
of density because a drop in the recombination cooling rate is compensated by a drop in the 
number of neutral atoms able to absorb energy by being photoionized. 

From Figure 1, we can already begin to guess how the impact of photoionization on a 
collapsing perturbation might depend on the parameters of the perturbation and of the UV 
background. Above T ~ lO^K, photoionization does not infiuence cooling curves at all, so 
we expect no impact on the formation of objects with virial temperature above this value 
(corresponding to a circular velocity Vdrc ~ 160kms~*, see equation [13] below). For a softer 
background spectrum, e.g. a = 5 instead of a = 1, the cooling bump at lO^K from singly ionized 
helium returns at moderate density, so the temperature above which cooling rates are minimally 
affected is lower. Furthermore, if (as we shall eventually argue) the heating of low temperature gas 
is more important than the reduction of cooling rates, then the relevant temperature to consider is 
Teg, which is considerably lower than the temperature at which cooling rates are unaffected. Again 
one expects a harder spectrum to have a greater impact. Finally, since equilibrium temperatures 
and reductions in cooling rates are smaller when the gas density is high, we expect photoionization 
to have less effect on objects that collapse at higher redshifts. 




(3) 
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2.2. Perturbations 



As initial conditions, we adopt the average density profile around a 2a peak in a Gaussian 
random density field of power spectrum Ak~^e~^ (Bardeen et al. 1986), as shown in Figure 2. 
The n = —2 power-law index is roughly the effective index of the cold dark matter power spectrum 
on galaxy mass scales. These initial conditions have two free parameters, A and Rf, which we can 
specify more physically by the filter mass Mj, defined as the mass contained within a sphere of 
radius 2Rf where Rf is the filter radius, and the collapse redshift Zc, defined to be the redshift at 
which the r = 2Rf shell would collapse to r = in the absence of pressure. We cast our definitions 
in terms of 2Rf instead of Rf because it is the radius at which the perturbation profile has fallen 
to roughly half of its peak value. 

We will often characterize a perturbation by the circular velocity 



= V ^ 

that would be associated with the shell surrounding mass Mf if it collapsed and virialized at half 
its turnaround radius r^^. The circular velocity can be expressed in terms of Mf and the collapse 
redshift Zc- The turnaround radius is given by 

rta = n (^^'^i^ 5 (5) 

where Si and are the mean interior density contrast and the radius at some redshift Zi high 
enough that the perturbation is still in the linear regime. The collapse time tc is given by 

where ti is the time corresponding to the redshift Zi, and the mass Mf is given by 

Mf = ^vrrf/>,(l + Zif, (7) 

where pc is the critical density of the universe today. From equations (5)-(7), and using the 
relation 



ti V 1 + z, 
we get 

M, 



The circular velocity can therefore be written as 



128kms-if^^^,^l (^] . aOl 




Fig. 2. — Mean profile of a 2(T peak in a Gaussian density field with power spectrum P{k) = Ak~^e~' 
The solid line shows the density contrast; the dashed line shows the mean interior density contrast, 
adopt this profile as the initial density distribution for our collapse simulations. 
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It is natural to express results in terms of a circular velocity rather than a perturbation 
mass because we can directly associate a temperature with this circular velocity. Let us assume 
that the local density profile is p(r) = pi(r/ri)~^ and demand that the gas be isothermal and in 
hydrostatic equilibrium. In this case, the gravitational force on a shell of width dr enclosing a 
mass M is balanced by the differential pressure force, 



4:TTr^p(r)dr = \--^dr]4:TTr\ (11) 
/- \ dr J 

where the pressure p is given by 

p=^^. 12 

lim,p 

Here, kB is the Boltzman constant, p, T and p are the pressure, temperature, and molecular 
weight of the gas, and is the proton mass. With Vdrc = \fGMjr^ equations (11) and (12) yield 

T = ( vL^ = 5.95 x10'k(A]( \ (13) 



— t;-. = t>.yt> X iU K — 

\2kB J VO.e; Vl28kms 

This result can also be obtained by equating the kinetic energy of a mass M moving at the rms 
velocity v^ms to the internal energy of a mass M at temperature T, 



For a singular isothermal sphere, the rms velocity is related to the circular velocity by 
Vrras = \/^T^'^^circ (Biuuey & Tremaiue 1987). Substituting into equation (14), we get 

and 

r=(f^)»:i. im 

the same result as given by equation (13). 



3. Results 



3.1. Physical picture 



To illustrate the effects of the UV background on the trajectories and temperature histories of 
fiuid shells, we performed runs for collapse redshift Zc = 2 and two values of the circular velocity, 
Vc,,,c = 35kms-i (Mj = 0.02 x lO^^M®) and Vc^rc = 128 km s'^ (Mj = lO^^M®). In each case 
we ran the code once without a UV background { J-21 = 0) and once with a redshift-independent 
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Fig. 3. — Fluid sliell trajectories in liigli mass collapses (Mf = 1 x 10^^ Mq, = 2, , v^rc = 128 km s""'^) 
with (dashed lines) and without (solid lines) a UV background. The time and radii are scaled with respect 
to the collapse time tg and the filter radius Rf. 




Fig. 4. — Time evolution of the temperature of a typical fiuid shell in the high mass collapses of Fig. 3, 
with (dashed lines) and without (solid lines) a UV background. 
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UV background of amplitude J -21 = 1 and spectral slope a = 1. The starting redshift for all 
simulations in this paper is Zi = 36. 

Figure 3 shows the trajectories of six widely spaced fluid shells from the high mass runs. 
The time and radii are scaled with respect to the collapse time tc and the filter radius Rj. The 
solid lines represent trajectories obtained with no UV background, while the dashed lines are the 
trajectories for J-21 = 1. There is little difference between the solid and the dashed trajectories. 
Figure 4 shows the time history of the temperature of a typical fiuid shell in the absence (solid line) 
and presence (dashed line) of the UV background. When the UV background is present, the gas is 
preheated to a temperature of about 25, 000 K. Since this temperature is an order of magnitude 
smaller than the temperature T ~ 6 X 10^ K associated with a circular velocity of 128 kms~^, this 
preheating does not affect the collapse of the gas. Therefore, as illustrated in Figures 3 and 4, the 
UV background has essentially no effect on high mass collapses, as we speculated in §2.1 on the 
basis of the cooling curves. The physics underlying the trajectories and temperature histories in 
the no-ionization case is discussed in Paper I. 

In Figure 5 we show fiuid shell trajectories from the low mass collapses. Once again, the solid 
lines are trajectories for no UV background, and the dashed lines are trajectories for J-21 = 1, 
a = 1. This time the trajectories are very different with and without the UV background. 
Figure 6 shows the temperature history of a typical fiuid shell in the absence (solid line) and 
presence (dashed line) of the UV background. When the UV background is present, the gas is 
again preheated to a temperature of about 25, 000 K. However, this temperature is now of the 
same order of magnitude as the temperature T ~ 45, 000 K associated with a circular velocity of 
35 kms~^. Pressure support from the warm gas delays turnaround substantially and slows the 
subsequent collapse, greatly reducing the amount of baryonic mass that is able to collapse, cool, 
and sink to the center in a given amount of time. The dynamics of the coUisionless dark matter 
shells (not shown) are almost identical with and without the UV background, since the baryons 
provide only a small fraction of the gravitating mass. 

Figures 5 and 6 show that photoionization can dramatically alter the evolution of a low mass 
perturbation. The fact that these effects appear before a gas shell shock heats and begins to cool, 
indeed before the shell even turns around, strongly suggests that the UV background exerts its 
infiuence by heating the pre-coUapse gas rather than by reducing cooling rates. The heavy dashed 
lines in Figures 5 and 6 confirm this view. They show trajectories and a temperature history from 
a run with J-21 = 1 in which we use photoionized abundances to compute cooling rates but set all 
heating terms to zero. The trajectories in this case are nearly identical to those obtained with no 
UV background at all. In this physically artificial scenario, cooling is possible at T < lO'^K because 
the photoionized gas emits recombination and Bremsstrahlung radiation and has no compensating 
heat input, so the gas cools well below lO'^K and shells collapse slightly faster than they do when 
J-21 = 0. 

Since the equilibrium temperature is lower for a softer UV spectrum, we expect the impact 
of the UV background to be smaller (see §2.1). The dotted lines in Figures 5 and 6 confirm this 
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Fig. 5. — Fluid shell trajectories in a low mass collapse (Mf = 0.02 x 10"'^"'^ Mq, Zg = 2, v^rc = 35kms~"'^). 
The solid lines are obtained in the absence of a UV background (J-21 = 0). The dashed lines are obtained 
with J_2i = 1 and a = 1. The dotted lines are obtained with J-21 = 1 and a = 5. The thick dashed lines 
are obtained with J-21 = 1 and a = I, but neglecting heating terms. The time and radii are scaled with 
respect to the collapse time tg and the filter radius Rf. 






12 3 

Fig. 6. — Temperature histories of a typical fiuid shell in the collapses of Fig. 5, with no UV background 
(solid line), J-21 = 1 and a = 1 (dashed line), J-21 and a = 5 (dotted line), and J-21 = 1, a = 1 with 
heating terms neglected (thick dashed lines). 
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expectation. They show results for the same perturbation (vdrc = 35kms~^, Zc = 2) and UV 
background amplitude ( J-21 = 1) but with spectral slope a = 5 instead of a = 1. Pressure 
support again delays turnaround and collapse, but the pre-heat temperature is lower and the 
resulting shell trajectories are intermediate between those obtained for no UV background and 
those obtained for a = 1. 

3.2. Suppression of low mass galaxy formation 

To study the implications of photoionization for galaxy formation, we have performed several 
series of collapse calculations. Within each series we vary the value of the filter mass Mj, and from 
one series to another we vary other parameters such as the collapse redshift and the amplitude, 
spectral slope, and time history of the UV background. In all of these calculations, we use ^l^ = 0.9 
and 1)5 = 0.1. In a hierarchical model of structure formation, a dark matter halo typically merges 
with another dark halo of comparable size in a time of order its own collapse time tc (Lacey & 
Cole 1993). Gas that has not cooled and condensed by the time a merger occurs will typically 
be shock heated to the virial temperature of the new, merged halo. As our primary diagnostic, 
therefore, we show the mass Mc of gas that cools and collapses to r = by time t = 2tc, and 3tc 
as a function of the circular velocity Vdrc- Here, tc = ^o/(l + ^cf'^'^ ^ where is the present age of 
the universe and is the collapse redshift as defined in §2.2. We scale to the mass Mdp = 0) 
of gas that would cool and collapse to r = in the same time if radiative cooling were perfectly 
efficient (i.e., pressure equal to zero). 

Figure 7 illustrates the effects of the UV background amplitude J-2i- It shows the ratio 
Mc/Mc(p = 0) for runs with no ionizing background (crosses), with J-21 = 1 (circles), and with 
J-21 = 0.3 (asterisks). In all cases Zc = 2, and in the last two cases the spectral slope is a = 1. 
As discussed already in §3.1, the photoionizing background inhibits the formation of low mass 
galaxies but does not change Mc in high mass collapses. Furthermore, there is little difference 
between the results obtained with J-21 = 0.3 and those obtained with J-21 = 1- In both cases, we 
find that baryon collapse is completely suppressed below a cutoff circular velocity Vdrc ~ 30 km s~^ 
at t = 3tc. At t = 2tc, Mc is reduced by about 50% for Vdrc ~ 50kms~^. 

To assess the effects of the shape of the UV spectrum, we compare runs with a = 1, a hard 
spectrum that might be expected if the UV background is produced mostly by quasars, and a = 5, 
a softer spectrum more like that expected from star-forming galaxies. In both cases we use an 
amplitude J-21 = 1 and a collapse redshift Zc=2. Figure 8 plots the results for no UV background 
(crosses), a = 1 (circles), and a = 5 (asterisks). The UV background inhibits the formation 
of low mass galaxies regardless of the value of a. However, the effect is greater for a = 1, as 
speculated in §§2.1 and 3.1. The cutoff circular velocity at t = 3tc is Vdrc ~ 30kms~^ if ci = 1 but 
Vdrc ~ 22kms~^ if Q! = 5. The collapsed mass at 2tc is reduced by about 50% at Vdrc ~ 54kms~^ 
for a = 1 and at Vdrc ~ 41 km s~^ for a = 5. 
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Fig. 7. — Mass Mc of gas that cools and collapses by time t = tc,2tc,3tc, as a function of the circular 
velocity Vrjrc in kms""'^. Here tc = to/{l + Zr)^/^, where to is the present age of the universe, and we have 
taken Zg = 2. We have normalized Mc to the mass Mc(p = 0) of gas that would cool and collapse in the 
same time in the absence of pressure. The crosses are obtained in the absence of a UV background; the open 
circles are obtained in the presence of the photoionizing flux defined by equations (1) and (2), with J_2i = 1 
and a = 1; the asterisks are obtained for J_2i = 0.3 and a = 1. 
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Fig. 8. — Like Fig. 7, but illustrating the dependence on the UV background spectral index a. Crosses 
show results for no UV background; open circles show results for a UV background with J_2i = 1 and a = I, 
and asterisks for J_2i = 1 and a = 5. 
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Fig. 9. — Like Fig. 7, but illustrating the dependence on the UV background time history. Crosses show 
results for no UV background, open circles show results for a constant photoionizing flux with J_2i = 1 and 
a = 1, and asterisks show results obtained with the time-dependent UV flux defined by equation (17). 
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The assumption of a constant UV background intensity from 2; = 36 to 2; = is unrealistic. 
To examine the effects of the UV background history, we compare the constant background 
( J-21 = 1) to a background that varies in roughly the way that might be expected if quasars are 
the dominant sources of the background: 



In both cases we adopt a spectral slope a = 1 and a collapse redshift Zc=2. The results for the 
two models, plotted in Figure 9 together with the no background case, are not very different. 
Suppression by the UV background is slightly weaker when we use the history of equation (17) 
instead of a constant J-21, since the gas in the former case is preheated by photoionization 
for a shorter time before turnaround. The effect is small, however, because the photoionizing 
background is turned on at a relatively high redshift, z = 6, and the amount of gas that turns 
around by this time is very small compared to the amount of gas that turns around by the collapse 
redshift Zc = 2. In both cases, the lower cutoff at t = 3tc is Vdrc ~ 30kms~^, and the collapsed 
mass at 2tc is reduced by about 50% at Vdrc ~ 50kms~^. 

In the case described above, the gas is assumed to be cold before the UV background turns 
on at 2; = 6, and the gas temperature thereafter evolves according to the heating and cooling rates 
predicted for ionization equilibrium. Miralda-Escude & Rees (1994) point out that the process of 
reionization can heat gas more rapidly than this if it occurs fast enough. To examine the possible 
impact of such rapid heating, we compare the last series of runs to a series with the same UV 
background history (equation [17]) but with the gas temperature set instantaneously to an initial 
temperature Ti = 50,000K at 2; = 6. The results are shown in Figure 10. The value of the gas 
temperature at 2; = 6 has no effect on the amount of gas that collapses by 2; = 2;c and later. The 
equilibrium temperature at 2; = 6 is not much below 50,000K, and the gas heats quickly to this 
temperature even from a cold start. 

Figure 11 illustrates the effect of the collapse redshift, comparing runs with Zc = 2 and 
= 5. In both cases J-21 = 1 and a = 1. The effect of the background is slightly weaker for 
the higher collapse redshift. The cutoff at t = 3tc drops from Vdrc ~ 30kms~^ for 2;c = 2 to 
Vcirc ~ 25kms~^ for Zc = 5, and the circular velocity for which Mc at t = 2tc is reduced by 50% is 
Vcirc ~ 54kms~^ for 2;c = 2 and Vdrc ~ 46kms~^ for Zc = 5. As predicted in §2.1, perturbations 
that collapse earlier are less affected by photoionization because their gas densities are higher and 
their equilibrium temperatures are correspondingly lower, but the difference is not striking. The 
difference would have appeared much larger if we had used mass Mf instead of circular velocity 
Vcirc as the independent variable in Figure 11, but as argued in §2, circular velcity is the more 
natural parameter for this problem because of its direct connection to temperature. 
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Fig. 10. — Like Fig. 7, but illustrating the effect of the gas temperature when the time-dependent UV 
background flux is turned on. Crosses show results for no UV background. Open circles and asterisks show 
results obtained with the time-dependent UV flux defined by equation (17), with an initial temperature at 
z = 6 of Tj- = and Ti = 50, OOOK, respectively. 
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Fig. 11. — Like Fig. 7, but illustrating the dependence on the collapse redshift Zg. Crosses and triangles 
show results for no UV background with Zg = 2 and = 5, respectively. Open circles (zg = 2) and asterisks 
(zc = 5) show results for a UV background with J_2i = 1 and a = 1. In all cases tc = to/{l + Zr)^/^, where 
to is the present age of the universe, and Vdrc is defined by equation (10). 
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4. Discussion 

Table 1 summarizes the results illustrated in Figures 7-11. For each series of runs, we list the 
circular velocity at which the cooled mass is reduced by 50% compared to its value in the absence 
of a UV background, and the cutoff circular velocity at which the baryon collapse is completely 
inhibited by the UV background, always at time t = 2tc. The main result is a reduction of about 
50% in the cooled mass at Vdrc ~ 50kms~^, and complete suppression at Vdrc ~ 25 — 30kms~^. 
As the table shows, this result is remarkably insensitive to the UV background parameters or 
collapse redshift. The most significant parameter is the spectral slope a, which determines the 
equilibrium temperature of the gas (see Figure 1). Because of the robustness of this result, we 
expect that it is also insensitive to some of the parameters that we have not varied, such as i7 and 
the shape of the initial density perturbation. 
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Table 1: Summary of the results at time t = 2tc for different values of the parameters: the collapse redshift 
Zc and the amplitude J-21, spectral slope a, and time history of the UV background. We list the circular 
velocity Vdrc at which the cooled mass is reduced by 50% compared to its value in the absence of a UV 
background, and the cutoff circular velocity at which the baryon collapse is completely inhibited by the UV 
background. 

Our study confirms Efstathiou's (1992) suggestion that photoionization can suppress the 
formation of low mass galaxies. The simulations yield quantitative estimates of the effect that 
are difficult to obtain in a robust way from Efstathiou's analytic approach. We are also able to 
clarify the physics of the effect: even though photoionization alters the cooling rates, it is clear 
from Figure 5 that it is the heating of the gas that matters. If the gas is able to overcome pressure 
support and collapse, it cools rapidly despite lowered cooling rates. Figure 5 also shows why 
the largest affected masses are somewhat higher than one might naively guess by equating the 
perturbation's virial temperature to the gas equilibrium temperature. Pressure support delays the 
turnaround of a gas shell, and the density at turnaround is lower than the density at virialization, 
so the characteristic dynamical temperature associated with the perturbation (which scales as 
T ~ \J M I p^l-^) is somewhat higher. 
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Our results are consistent with the three-dimensional, SPH simulations of Weinberg, Hernquist 
& Katz (1995), who find no suppression of galaxy formation at 2; = 2 down to their resolution 
limit, which corresponds to a circular velocity Vdrc ~ 100 kms~^. Our results imply that their 
minimum resolved mass is still substantially above the mass where photoionization infiuences 
galaxy formation. Still more reassuring, our results are consistent with the high-resolution SPH 
simulations of individual collapses by Quinn, Katz, & Efstathiou (1995). They find a drop in the 
cooled baryon mass of ~ 30% at Vdrc ~ 45kms~^, and complete suppression of baryon collapse 
at Vcirc ~ 25kms~^, values remarkably similar to those in Table 1. The two studies complement 
each other nicely. Quinn et al.'s results indicate that our conclusions are not sensitive to the 
idealizations of spherical symmetry and smooth initial conditions. Our results indicate that the 
suppression found by Quinn et al. is not sensitive to the numerical resolution (inevitably lower 
in 3-d simulations) or to their assumed parameters of the UV background. In particular, since 
the effects of the UV background depend mainly on the equilibrium temperature, the strong 
background intensity ( J-21 = 10) assumed by Quinn et al. should not make much difference to 
the scale of suppression. Vedel et al. (1994) find that the internal structure of Milky Way-size 
galaxies can be affected by the presence of a UV background. Although they attribute this effect 
to the lower cooling rates, it is more likely a result of the reduction in small scale substructure 
during the collapse because of gas heating. 

By suppressing collapse into low mass objects, photoionization can solve one problem that 
has faced hierarchical models of structure formation, a tendency to cool nearly all gas into small 
systems at very high redshifts. If first generation stars ionize the surrounding intergalactic medium, 
they will delay the formation of the next generation of objects until perturbations larger than the 
Jeans mass at T ~ few X lO'^K are able to collapse. It is less clear that photoionization can solve 
the problems of hierarchical models in reproducing the observed galaxy luminosity function, since 
Press-Schechter (1974) type calculations indicate a rise above the estimated luminosity function 
already at Vdrc ~ lOOkms"^ (White & Frenk 1991). However, photoionization will certainly help 
to reduce the discrepancy, and estimates of the slope of the local luminosity function may go up 
as low surface brightness galaxies are taken into account more accurately. On the theoretical 
side, a natural next step is to incorporate the results of collapse calculations like those presented 
here into Monte Carlo models of galaxy formation history (Kauffmann et al. 1993; Cole et al. 
1994; Heyl et al. 1995). The inclusion of realistic photoionization and cooling calculations will 
solidify the physics underlying the Monte Carlo models, and the models will provide a larger 
cosmological framework within which to place the results of individual collapse simulations. The 
combination will allow better comparisons between the galaxy population predicted from theory 
and the galaxies in the local and high-redshift universe. 
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